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Inelastic spin relaxation and spin splitting Es in lateral quantum dots are studied in the regime 
of strong in-plane magnetic field. Due to both g-factor energy dependence and spin-orbit coupling 
Es demonstrates a substantial non-linear magnetic field dependence similar to that observed by 
R.Hanson et al [Phys. Rev. Lett. 91, 196802 (2003)]. It also varies with the in-plane orientation 
of magnetic field due to crystalline anisotropy of spin-orbit coupling. Spin relaxation rate is also 
anisotropic, the anisotropy increasing with the field. When the magnetic length is less than the 
'thickness' of GaAs dot, the relaxation can be order of magnitude faster for B|| [100] than for B|| [110]. 
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Proposals to use electronic spin in quantum dots for 
quantum information processing have fuelled extensive 
studies of spin-orbit (SO) coupling in heterostructures as 
means to manipulate the electron spini and as a source of 
spin relaxation. Recent theories^iii and experiments^iSii 
suggest that spin relaxation in quantum dots is strongly 
suppressed by electron confinement but may be sped-up 
by a magnetic field, in particular, by the field parallel to 
the plane of the lateral stucture. 

It is customary to assume that an in-plane magnetic 
field couples only to spin of the electron. In this ap- 
proximation one can describe spin relaxation in terms 
of effective two-dimensional (2D) SO coupling2ii4. This 
aproach can be justified provided that Xb > A^, where 
Xb ~ \/ ch/eB is magnetic length and Xz is the extent 
of the subband wave function across the 2D plane. In 
the oposite limit that corresponds to a strong magnetic 
field, Xb "C Xz, subbands in a heterostructure trans- 
form into bulk Landau levels (magneto-subbands) thus 
changing parameters of the effective 2D motion^. This 
effect has been observed in optical and FIR spectroscopy 
of low-density GaAs/AlGaAs heterostructures^, resonant 
tunnelling in double-barrier devices^, and in quantum 
transport characteristics of lateral dotsii. 

In this Letter, we propose a theory of the spin relax- 
ation of electrons in lateral dots in a strong in-plane mag- 
netic field. The field effect on the orbital electron motion 
transforms states in low-density heterostructures into 
magneto-subbands. We take into acount this crossover 
as well as the Dresselhaus-type spin-orbit coupling in 
GaAalS. We show that at high fields both the inelas- 
tic spin-flip time Ti at low temperatures kT <^C £s and 
the electron spin splitting Eg depend on the magnetic field 
orientation with respect to crystallographic axes, which 
can be used to distinguish the SO coupling-induced ef- 
fects from those caused by a hyperfine interaction with 
nuclei. We also present analytical description of £s and 
Ti dependences on both the magnitude and direction of 
the magnetic field. 

The Hamiltonian of electrons in a dot made of a lateral 
GaAs/AlGaAs structure grown in direction 1^ = [001] 



can be written as 



3D 



pI+Px+ Py 



2m 



-3 



(1) 



kij—x,y,z 



In Eq. ^ we use two systems of in-plane coordi- 
nates. Axes X and y (used in the SO coupling term 
Hso) are determined by cristallographic directions [100] 
and [010], respectively. Axis X is directed along the in- 
plane field B = IxB, with 1^ = {lx,ly,0) and Y along 
ly — {—ly,lx,0)- In Eq. Q the kinematic, Pa = —ihda 
and canonical, Pa momenta are written in the coor- 
dinate system X and Y. We use the Landau gauge 
A = — (z — a)_Bly, so that Pz — Pz, Px — Px and 
Py — Py — ^B{z — a), with a to be specified later. 

In the spin part of iJaD, is Bohr magneton, g-factor in 
GaAsi^ is g ~ — 0.44-|-p^g', e'^*' is the antisymmetric ten- 
sor, and SO coupling constant 7 according to Refsii^ii^*i& 
3 

is 7 = (26 ± 6) eVA . Hso in Eq. ^ is written in the 
form which guarantees that it is Hermitian, despite the 
non-commutativity of operators px and py with p^pi. 

The dot is formed by a potential profile V = Vz{z) + 
^m'd'^{X^ + y^), which is stronger in the heterostruc- 
ture growth direction 1^ than within the XY plane. We 
consider two particular cases: triangular well Vz = Fz 
and parabolic well Vz{z) — ^muj'^z^. The wave functions 
\n,px,PY) of electrons in the n-th magneto-subband and 
their 2D dispersion, £n{PY) + '^Px determined by 



Hz 



^pl + Vz{z) + \m^l 



^^^Py] 



(2) 



whereas the parameter a is chosen in such a way that the 
lowest subband dispersion £o(py) has minimum at py = 



0. This defines z = 



and 



For a triangular well, we describe magneto-subbands us- 
ing the function of a parabolic cylinder^'* and evaluate 
a and my = rjm numerically. For a parabolic well, har- 
monic oscillator functions give a = and -q ~ \ + lo'I/ljj'^ . 

Since electron confinement across the plane is much 
stronger than in lateral directions, px,PY ^ Pz, we 
substitute py — py — -Bz = py — mUJcZ, Px = 



2 



IxPx - lyPY + lymuocz and py = lypx + IxPy - IxmcUcZ 
into -ffsD, expand it in powers of kinematic momenta 
px and Py (up to quadratic terms) and derive the ef- 
fective 2D Hamiltonian H2b{px,Py,'^)- Iu particular, 
when analysing SO coupling, we expand Hgo up to linear 
order in py and px , 
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where 



(3) 



Ox 



PzZPz 



\(pI 




1 Pxcrx 
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ZPzZ 



Pz Py<Jy 



Ixly 



X4 



Pxcty 

h 

Py<Jx 



{[ll-ll]pY + 2lJyPx) 



ZPz 



PzZ 



In both TJgQ and Hl^ the last term does not contribute 
to the effective 2D Hamiltonian: for magneto-subbands 
determined by in Eq. ^ {Q,px,Py\zPzz\0,pxtPy) = 
and {Q,px,Py\zPz+Pzz\^,Px,Py) = 0. 

The first term in //g^ yields an anisotropic addition 
to the 2D electron spin splitting linear in 7. The sec- 
ond term slightly turns the spin quantization axis off the 
magnetic field H = \xB. It can be neglected as long as 
we restrict ourselves by the lowest order in 7. Thus, 



Es = gfiB - IxlyjX^ As, 

g « -0.44+ (Ob^lO)g', A, - (hymF) 
A3(0|53|0) ^Xl{0\pz~ZPz\0) 



(4) 



1/3 



As 



Xlh^ 



The anisotropy of spin splitting is crucially sensitive to 
the inversion asymmetry of the confinement potential 
Vz, thus it is a peculiarity of heterostructures. The 
anisotropy effect in Eq.l^J is maximal in a field oriented 
along crystallographic directions [110] or [110]. The field 
dependence of the anisotropic part of spin splitting is 
characterised by the parameter As. In a weak mag- 
netic field, ujcfi < £1 — £o: perturbation theory analy- 
sis gives As « 2A6^^^^^lUi. = ZA2lijJi/ {ei — £0) leading 
to the anisotropy of linear g- factor. The field depen- 
dence As{B) at high fields is shown in Fig. 1(a). For 
GaAs/AlGaAs heterostructure with Xz ^ lOOA and 

a 3 

7 ~ (26±6)eVA , Eq.lO predicts that the spin split- 
ting £s is modulated by about 10% for different orienta- 
tions of the magnetic field. Ss also includes an isotropic 
non-linear B-dependent part due to the g-factor depen- 
dence on the electron momentum, (0|p^|0) ~ Beh/2c, 
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(a) 



As 



my /rn 




XieB/Tic 



FIG. 1: Magnetic field dependence of parameters (a) As, rj = 
my /m and (b) Qx,y, ^x,y for potential well 14 = Fz. 



thus g^iB K. -O.MfiB + ^g'^B'^. A non-linear £s(B)- 
dependence similar to that described by Eq. |0J was 
reported in Refi^ where measurements have been made 
in a field applied along [110] axisSi. 
In the effective 2D Hamiltonian, 



H2U — 



jpx - ax) jPY - o-y) 
2m 



2r]m 
-Ox, 



(5) 
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ax = -fiXso^ {[ll - ll] Oxcrx + IxIvOycty) , 

SO coupling arises from the first three terms in Hg^. It 
manifests itself via non- Abelian gauge fields ax , ay which 
display an anisotropy linked both to the direction of ex- 
ternal magnetic high field B = Ijs^i? and crystalline axes, 

Ox = fiXl(K - 0, Qy - -hXl{2K + 0, 
Tx = ~hXl{2K - 3C)ry, Ty ^ -hXlnrj, 



Fig. 1(b) shows how O and F depend on the mag- 
netic field^® for a triangular well Vz = Fz [with Xz = 



3 



1 /3 

{ji^ /mF) ]. At high fields these dependences are simi- 
lar to what we found for a parabolic well Vz{z) = ^muj^z^ 
with Az = yj2hlmLo, 

Ox = -■:Qy = --5g;rx = g^-Sg-Sy = -g^ 
g g 



where g ~ ^yT+'ujfJuj^ « XleB/2hc at lUc S> uj- This 
similarity implies that at high field we can approximate 
&, r, and TOy in hctcrostructures by their values obtained 
for a parabolic well with the same Xz ■ 

Lateral orbital states described by H2D have 
the spectrum Emm' = (M + + [M' + 

^)h'dr]~'^/'^ . The lowest level wave function^® is |0) — 

{TrXXyr^^^ e-^'^^'e-'^'^^'^ifaiz), where A = ^2hlmd 
and Ay = 77^^/^A. Here, -d and r]~^^^'d are the fi'cquen- 
cies of electron harmonic oscillations along the X and Y 
axes, respectively, and the dot states |n) = |M, M' ) are 
characterized by quantum numbers M and AI' . 

The rate of the phonon-assisted spin fiip |0, +) 
0,—) in the lowest order in both the e-ph interaction 
and SO coupling is 



A = 



L^dq 
0\W\n){n\hl\0] 



27r 



\AfS{es~hsq), (6) 
. , (0|/iro|n)(n|T4-|0)- 



Eq — En+ Es 



Eq — En — Es 



Here, W{r, q) — w e^'^^'/L'^/^ is the phonon field with 
being the normalisation volume for phonons. We choose 



w 



q. 



where 9s = 7— • 
ns 



to take into account both piezoelectric and deforma- 
tion (S) phonon potential. 

Operators /i"o(P; i") can be obtained using Eq. 



hence it is sufficient to analyse A using only the linear in 
momentum SO coupling in i?2D- 

In a parabolic dot operators px and py couple the 
state |0) = |0,0) only to states |0, 1), |1,0), and 

(0,0|e^'i"-|l,0) = (l,0|e''i--|0,0)-^gxAA, 
(0,0|e''i"-|0,l) = (0,l|e''i--|0,0) - ^gyAyA, 

where f{qz) = J dze^'^'^\(pQ{z)\'^ . As a result, 

_ wes AA^ [ IJyQyqx ^ [^l - ^l] Tygy \ ^^^^ 



M 2Xso 1 1 - {e,/My 1 - {bs/MY 



V 



The angular distribution of the phonon emission is de- 
termined by the form-factor A(q) in Eq.®. Depending 
on the magnetic field, the emitted phonon wavelength 
qs — Ss/hs may fall into one of the following regimes: 



A) q, < A-i; B) A-i <q,< Ay^ C) A^^ < 



In regime A, the phonon wavelength excedes all dimen- 
sions of quantum dot. As a result, Awl and phonons 
are emitted isotropically. In regime B, most of phonons 
are emitted perpendicularly to the magnetic field direc- 
tion, since the phonon wavelength is shorter than the lat- 
eral dot size A in the direction of external field. Accord- 
ingly, e"}-^ w 1, e-(«^^/2)' « 1, thus / « 1 and jAp « 
g-(9xA/2) Finally, in the high- field regime C phonons 
are emitted across the heterostructure. Using the simi- 
larity between magneto-subband states and bulk Lan- 



dau levels, we approximate / « exp 



and 



|A|^ 



exp 



- {qxX/2f - {qyXy/2f - i feAs)' 



The rate can be evaluatedSSi using Eqs. ©-(CDIl, 



a ■ /iso(p,r) 



Pxax pyay 



m 



my 



(7) 



As long as the orbital part of the Hamiltonian iJ2D 
remains T-invariant, the orbital eigenstates are real, 
and (0|e*i'"|n) = (n|e*i""|0). Moreover, {n%o\Q) = 
— (O|/iso|n) because spin operator a changes sign under 
the t —t transformation whereas the product S-hso re- 
mains the same (as a spin-orbit part of T-invariant i?2D)- 
Consequently, two terms in the amplitude of the phonon- 
emission-assisted spin-flip process in Eq. © cancelSi in 
the limit £s — > 0, and the transition amplitude reads 



A = 2wes ^ 



(0|/i^ \M,M'){M, M'|e'i''|0) 



M,M'>1 {h-dM + M'M^/rn/rn^ 



(8) 



Being generic for any T-invariant Hamiltoniar^ such a 
cancellation should take place in all orders in px and py. 



M 



Qx 



(11) 



(aljy) 



\ [1 - Tiel/h?d-^Y [1 - el/hH-^Y J 

where w1 — {/3'^/qs) -f- S^^s, while Q and a are specific 
for each particular regime (A-C), 



?B = ^ [Xqsf , 

Qc = (A/Ay)3e-5(9=^«)', ac = QyXy/TyX. 



UA = Oy/Fy; 
Qb = (A(?s)^ , QfB = 2ey/9sAry; 



The factor in curly brackets in Eg. 1)11(1 determines the 
relaxation rate dependence on the magnetic field orienta- 
tion. If the field is so weak that Xz < Xb, then ^ « 2, 
A ~ Ay, a — 2, and T^^ turns out to be isotropic. 
The anisotropy develops when Xb ^ Xz [i.e., 77 > 1, 
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A > Ay and Qy/^y < 2] and increases with the field. 
According to Eq. l(TT|l. at high fields where A^ ^ A2 
but £s < M spin relaxes faster in a magnetic field ori- 
ented along [100] or [010] and slower when B is paral- 
lel to [110] or [110], which was the field oricntationSi 
in the experiment in Ref.^. In the field range where 
As ^ Az and Es < M/rj^^^, spin- flip rate for those two 
orientations has power-law different field dependences^, 
Tf 1 (B||[100]) oc 517/2 and Tf ^ (B||[110]) oc B^^. 

The anisotropy in Ti is strongly enhanced in the vicin- 
ity of crossing of the level |0, with |0, 1, — ) or |1, 0, — ), 
though the divergence ofT^^ in Eq. l(TT|l at Sg = M/rj^^'^ 
and £s = h-d is an artifact of the lowest-order pertirbation 
theory analysis and it is prevented by level anti-crossing 
due to SO coupling22,. For B||[100] or B||[010], spin re- 
laxation is resonantly sped- up when Eg — hd/'q^/'^. For 
B||[110] and B||[110] the rate acquires a maximum 
at a higher field where Es — hd. For samples used in Refi^ 
M « ImeV, thus the crossing of |0, and |1, 0, — ) levels 
was beyond the experimental field range. However, for 
Aj ~ lOnm and M ~ ImeV the crossing of levels |0, 4-) 
and |0, 1, — ) should enhance spin relaxation at B around 
15 ^ 20T if B II [100] or B II [010]. The formula in Eq.ini) 



is not exact when K'djri^l'^ < Eg < Wcg, nevertheless, 
in that field range the anisotropic behavior of T^^ (B) 
persists, since the spin-flip for B||[100] or B||[010] is en- 
hanced due to the openning of additional relaxation chan- 
nel |0,+) ^ |0,1,-). 

To conclude, we studied the effects of the spin-orbital 
coupling on the spin splitting £s and inelatic spin relax- 
ation rate T^^ in lateral quantum dots at low tempera- 
tures kT <C £s- We found that Eg demonstrates a sizeable 
non-linearity and anisotropy in its fleld-dependence, Eq. 

The anisotropy in the spin relaxation, Eq. ifTT)) is 
predicted to be even stronger: if the magnetic field B 
is high, Xb ^ Az, the relaxation can be order of mag- 
nitude faster for B||[100] than for B||[110]. The latter 
feature of the spin relaxation due to SO coupling can be 
used to distinguish it from the spin relaxation involving 
hyperfinc interaction with nuclei. 

We thank L.Kouwenhoven, D.Loss and C.Marcus for 
useful discussions. This work was supported by the 
EPSRC-Lancaster Portfoho Partnership, INTAS 03-51- 
6453, ARO/ARDA (DAAD19-02-1-0039) and DARPA 
under QuIST programme. 



1 H.-A. Engel et al, Phys. Rev. Lett. 93, 106804 (2004); G. 
Burkard and D. Loss, Phys. Rev. Lett. 88, 047903 (2002) 

2 A.Khaetskii and Y.Nazarov, Phys. Rev. B 61, 12639 

(2000) ; Phys. Rev. B 64, 125316 (2001) 

3 I.L. Aleiner and V.I. Fal'ko, Phys. Rev. Lett. 87, 256801 

(2001) ; J.-H. Cremers, P.W. Brouwer, and V.I. Fal'ko, 
Phys. Rev. B 68, 125329 (2003) 

* V.Golovach, A.Khaetskii, D.Loss, Phys. Rev. Lett. 93, 
016601 (2004) 

^ R.Hanson et al, Phys. Rev. Lett. 91, 196802 

(2003); J.M.Elzerman et al, Nature 430, 431 (2004); 

J.M.EIzerman et al. Appl. Phys. Lett. 84, 4617 (2004); 

M.Kroutvar et al. Nature 432, 81 (2004) 
^ T.Fujisawa et al. Nature 419, 278 (2002) 

D.M. Zumbiihl et al, Phys. Rev. Lett. 89, 276803 (2002) 

* In heterostructures with 2 x lO^ielectrons per cm^, the 
reconstruction of the 2D electron spectrum has been 
observed— in the field range 5 7T, whereas at the elec- 
tron density 10^°cm~^ [i.e., one electron in a lOOOA-size 
quantum dot] magneto-subbands would already form at 
B ~ 2 3T. 

^ I. Kukushkin et al, JETP Lett. 51, 436 (1990); B.Meurer, 
D.Heitmann, K.Ploog, Phys. Rev. B 48, 11488 (1993) 

C. Kutter et al, Phys. Rev. B 45, 8749 (1992) 

D. M. Zumbuhl et al, Phys. Rev. B 69, 121305 (2004) 
F.Malcher, G.Lommer, U.Rossler, Phys. Rev. Lett. 60, 729 
(1988); R. Eppenga and M. Shuurmans, Phys. Rev. B 37, 
10923 (1988) 

M.J.Snelling et al, Phys Rev. B 44, 11345 (1991); 
R.M.Hannak et al. Solid State Commun. 93, 3132 (1995) 
B.Jusserand et al, Phys. Rev. Lett. 69, 848 (1992); 
B.Jusserand et al, Phys. Rev. B 51, 4707 (1995); 
L.Wissinger et al, Phys. Rev. B 58, 15375-15377 (1998) 
J.B. Miller et al, Phys. Rev. Lett. 90, 076807 (2003); W. 



Knap et al, Phys. Rev. B 53, 3912-3924 (1996) 
In low-density structures the Rashba coupling 
[Y.Bychkov and E.Rashba, JETP Lett. 39, 78 (1984)] is 

■n .12 14 

msignmcan^— 

Px = IxPX — lyPV, Py = lyPX + IxpY , CFx = IxCTX — lyO^Y 

and ay — lyOx + Ix^y- 

jn,0,py) are functions of parabolic cylinder 
Z-^U {-E, V2 [5Agi - pYXsh-"^ + A|AJ=*] ) , where 
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